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We study the dynamical instability of a spherically symmetric anisotropic fluid which collapses 
adiabatically under the condition of vanishing expansion scalar. The Newtonian and post Newtonian 
regimes are considered in detail. It is shown that within those two approximations the adiabatic 
index Fi, measuring the fluid stiffness, does not play any role. Instead, the range of instability 
is determined by the anisotropy of the fluid pressures and the radial profile of the energy density, 
independently of its stiffness, in a way which is fully consistent with results previously obtained 
from the study on the Tolman mass. 
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I. INTRODUCTION 



The problem of stability is of the utmost relevance in 
the study of Newtonian and general relativistic models 
of self-gravitating objects. This becomes evident if wc 
recall that any static stellar model, in order to be of any 
use, has to be stable against fluctuations. Furthermore, 
such a problem is closely related to the one of structure 
formation, since different degrees of stability/instability 
will lead to different patterns of evolution in the collapse 
of self-gravitating objects. It is therefore not surprising 
that a great deal of work has been devoted to this issue 
since the pioneering paper by Chandrasekhar 

Extensions of Chandrasekhar's result to non-adiabatic 
fluids 3, __|, anisotropic fluids [_| and shearing viscous 
fluids 5j, have been carried out in the past. In all of 
these works the key variable is the adiabatic index T\, 
whose value defines the range of instability. Thus for a 
Newtonian perfect fluid, the system is unstable for Ti < 
4/3. In the above mentioned references it was shown how 
different physical aspects of the fluid affect the range of 
instability of the system. More recently Q the stability of 
anisotropic stars with quasi-local equation of state @, __| 
has been invetigated. 

This work represents another forward step in that di- 
rection. Our main goal here consists in studying the dy- 
namical instability of a spherically symmetric fluid dis- 
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tribution, under the assumption of vanishing expansion 
scalar. 

The main motivation to undertake such an endeavour 
is provided by the following argument: Highly energetic 
explosions in self-gravitating fluid distributions are com- 
mon events in relativistic astrophysics (see for example 
[__, [l_| and references therein). Accordingly, a relevant 
question related to this issue is 

• How the system evolves after the explosion? 

Now, in a recent series of paper it has been stablished 
that the expansion-free condition is particularly suitable 
for describing that kind of phenomena [lll - [l3j |. 

Indeed, since the expansion scalar describes the rate of 
change of small volumes of the fluid, it is intuitively clear 
that the evolution of an expansion-free spherically sym- 
metric fluid distribution is consistent with the existence 
of a vacuum cavity within the distribution which would 
be formed after a central explosion. This point is further 
discussed in this paper. Therefore potential applications 
of our results are expected for those astrophysical scenar- 
ios where a cavity within the fluid distribution is likely 
to be present (see for example Furthermore, as the 

fluid (under the expansion-free condition) reaches the 
centre entails a blowup of the shear scalar, whose conse- 
quences in the app earance of a naked singularity cannot 
be overlooked (15f . 

The first known model sa tisfy ing the expansion-free 
condition is due to Skripkin [1 61 ] for the particular case 
of a non-dissipative isotropic fluid, with constant energy 
density. In that work, Skripkin addressed the very inter- 
esting problem of the evolution of a spherically symmet- 
ric fluid distribution following a central explosion. As a 
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result of the conditions imposed by him a Minkowskian 
cavity should surround the centre of the fluid distribu- 
tion. 

Only recently this model has been again addressed. In 
["TO a general study on shearing with vanishing expansion 
scalar of spherical fluid evolution is presented, which 
includes pressure anisotropy and dissipation. While in 
[l2| it is shown that the Skripkin model is incompatible 
with Darmois' junction conditions [l7|, and that inho- 
mogeneous expansion-free dust models are deprived of 
physical interest since they imply negative energy den- 
sity distributions. Further analytical solutions describ- 
ing expansion-free evolution may be found in 13]. Cav- 
ity evolution under kinematical conditions other than 
expansion- free, have been considered in |18j . 

From the above reasons, it follows that most, phys- 
ically meaningful, expansion-free models should require 
anisotropy in the pressure and energy density inhomo- 
geneity. Thus we shall consider in this paper locally 
anisotropic fluids (further arguments to .ju stify such kind 
of fluid distributions may be found in [19l - |2l| and refer- 
ences therein). Also, even though it is already an estab- 
lished fact, that gravitational collapse is a highly dissi- 
pative process (see [22|, HH and references therein), we 
shall restrain here for simplicity to adiabatic evolution. 
The role played by different dissipative processes in the 
instability of a fluid distribution have been discussed in 
detail in [2|-|5j]. 

The fluid under consideration has two delimiting hy- 
persurfaces. The external one separating the fluid distri- 
bution from a vacuum Schwarzschild spacetime and the 
internal one, delimiting the cavity within which there is 
Minkowski spacetime. Thus, we have to consider junc- 
tion conditions on both hypersurfaces. 

It should be mentioned that for cavities with sizes of 
the order of 20 Mpc or smaller, the assumption of a spher- 
ically symmetric spacetime outside the cavity is quite 
reasonable, since the observed universe cannot be con- 
sidered homogeneous on scales less than 150-300 Mpc. 
However for larger cavities it should be more appropriate 
to consider their embedding in an expanding Lemaitre- 
Friedmann-Robertson- Walker spacetime (for the specific 
case of void modeling in expanding universes see (24 . [25[ 
and references therein). 

Next we shall present our perturbative scheme which is 
very similar to that employed in but now with the 

additional expansion-free condition. The study of the 
resulting equations shows that at Newtonian and post 
Newtonian approximations, the range of instability is in- 
dependent on Ti. This result is intuitively clear if we 
recall that, on the one hand the expansion- free condition 
implies, roughly speaking, that (at least close to New- 
tonian regime), the fluid evolves without being "com- 
pressed", and on the other that Ti somehow measures 
the (in) compressibility of the fluid. Therefore it is to be 
expected that not very far from Newtonian regime T\ 
does not appear in the discussion under expansion-free 
condition. Beyond post Newtonian approximation, intu- 



ition is no longer a reliable guide and T\ is expected to 
play a role in the discussion. 

The plan of the paper is as follows. In section II we 
give the energy-momentum tensor, the field equations 
and the junction conditions. Section III is devoted to 
explain the physical implications of the expansion-free 
evolution. The perturbation scheme is presented in sec- 
tion IV. In section V a brief discussion on Newtonian 
and post Newtonian approximations is presented. The 
dynamical expansion-free equation is obtained in section 
VI. Finally there is a concluding section. 



II. THE ENERGY-MOMENTUM TENSOR, THE 
FIELD EQUATIONS AND JUNCTION 
CONDITIONS 

We consider a spherically symmetric distribution of 
collapsing fluid, bounded by a spherical surface E( e ).The 
fluid is assumed to be locally anisotropic. Choosing co- 
moving coordinates inside Y\^ e \ the general interior met- 
ric can be written 

ds 2 _ = -A 2 dt 2 + B 2 dr 2 + R 2 (d9 2 + sin 2 0d<p 2 ), (1) 

where A, B and R are functions of t and r and are as- 
sumed positive. We number the coordinates x° = t, 
x 1 = r, x 2 = 9 and x 3 = <f>. Observe that A and B 
are dimensionless, whereas R has the same dimension as 
r. Two radii are determined for a collapsing spherical 
fluid distribution by the metric ((T]). The first is deter- 
mined by R(t, r) representing the radius as measured by 
its spherical surface, hence called its areal radius. The 
second is obtained out its radial integration J B(t,r)dr, 
hence called proper radius. These two radii, in Einstein's 
theory, in general need not to be equal unlike in Eu- 
clidean geometry (fluids with both radii equal are studied 
in HI). 

The matter energy-momentum T~p inside £( e ) has the 
form 

T «f> = 0* + P±)V a Vp + P±g a p + (P r - P±)x a xp, (2) 

where /x is the energy density, P r the radial pressure, 
P± the tangential pressure, V a the four-velocity of the 
fluid and \a a unit four-vector along the radial direction. 
These quantities satisfy 

V a V a = -l , X a X a = l , X a V a =0. (3) 

The four-acceleration a a and the expansion O of the fluid 
are given by 

a a = V aiP V f> , & = V a , a , (4) 
and its shear o~ a p by 

o-a/3 = V( a -i3) + a( a Vp) - ^Q(g a p + V a Vp). (5) 
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Since we assumed the metric ([I} comoving then 

V a = A~ 1 5q, x a =B- 1 6«. (6) 

From ^ with ((5]) we have the non zero component for 
the four-acceleration and its scalar, 



A' a (A' 
° 1 = T' a aa= \AB 



and for the expansion 
9 = 



- - 

A \ B 



' R 



(8) 



where the prime stands for r differentiation and the dot 
stands for differentiation with respect to t. With ([6]) we 
obtain for the shear ([5]) its non zero components 



and its scalar 



'-B 2 a, 022 = 033 sin 2 I 



C (Ja/3 = gCT 



f 



-R J a, 



where 
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(9) 
(10) 

(11) 



Einstein's field equations for the interior spacetime (fTJ) 
to £' e ' are given by 

(12) 



G aj 3 — K T aj3 , 



and its non zero components with (TTJ) , and © become 
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.(16) 

The component (Tl4")) can be rewritten with © and (fTTj) 

as 

1 



4- 



(17) 



The mass function m{t, r) introduced by Misner and 
Sharp [23] (see also [28j]) reads 



R o 23 R 

m = -R 23 2 * = 2 



R 



(18) 



(7) and the non trivial components of the Bianchi identities, 



T~ a0 = 0, from (I2J) yield 



;/3 



or using © 



1 



A+(M + ^)| + 2( Ai + P ± )| 



i? 



A + (m + Pr)^e + 2(p ± - P r )i = 0, 



and 



T~ a/3 v - 



B 



= 0, 

(19) 
(20) 

= 0. 
(21) 



Introducing the proper time derivative Dt given by 

1 d 



D T = 



Adt' 



(22) 



we can define the velocity U of the collapsing fluid (for 
another definition of velocity see section III) as the vari- 
ation of the areal radius with respect to proper time, i.e., 

U = D T R. (23) 

Using field equations and (I18p we may write 

/;/ = ^nR'R 2 . (24) 



Outside £( e ) we assume we have the Schwarzschild 
spacetime, i.e. 

ds 2 = - f 1 - ?^pj dv 2 - 2dpdv + p 2 {dd 2 + sin 2 9d<f> 2 ) 7 

(25) 

where M denotes the total mass, and v is the retarded 
time. 

The matching of the adiabatic fluid sphere to 
Schwarzschild spacetime, on the surface r — r S ( e ) = con- 
stant (or p — jo(w) S (e> in the coordinates of (1251) ) , re- 
quires the continuity of the first and second differential 
forms (Darmois conditions), implying 



rW / 2M 
Adt = dv I 1 



R = p(v) 



m(t,r) S =' M, 



(26) 

(27) 
(28) 



and 



where = means that both sides of the equation are eval- 
uated on £< e ). Comparing $Z9§ with (JT3J} and (JH]) one 
obtains 



Pr 



0. 



(30) 



Thus the matching of ([T]) and (I2~5l) on £( e ) produces 
and ([501) . 

As we mentioned in the introduction, the expansion- 
free models present an internal vacuum cavity (reasons 
for the formation of this cavity are discussed in the fol- 
lowing section.). If we call £W the boundary surface 
between the cavity and the fluid, then the matching of 
the Minkowski spacetime within the cavity to the fluid 
distribution, implies 



m(t,r) = 0, 



Pr 



(31) 



(32) 



III. ON THE PHYSICAL IMPLICATIONS OF 
THE EXPANSION-FREE EVOLUTION 



In section II we introduced the variable U which, as 
mentioned before, measures the variation of the areal ra- 
dius R per unit proper time. Another possible definition 
of "velocity" may be introduced, as the variation of the 
infinitesimal proper radial distance between two neigh- 
boring points (SI) per unit of proper time, i.e. Dt(SI). 
Then, it can be shown that (see [ll| for details) 



D T (Sl) 
51 



e), 



or, by using §8$ and (fTTj) . 

D T (Sl) 
51 



B 

AB' 



Then with ©, (HI]), ([23]) and ([34]) wc can write 



D T (6l) D T R _ D T (5l) 



SI 



R 



51 



and 
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D T (5l) 2D T R D T {61) 



51 



R 



51 



U 
R : 

2U 
~R' 



(33) 



(34) 



(35) 



(36) 



Thus the "circumferential" (or "areal") velocity U, is re- 
lated to the change of areal radius R of a layer of matter, 



whereas Dt{51), has also the meaning of "velocity", be- 
ing the relative velocity between neighboring layers of 
matter, and can be in general different from U . 

We shall see now that the condition = is associated 
to the existence of a cavity surrounding the centre of the 
fluid distribution. Indeed, if = then it follows from 
flSSJ) and flU that 



a= - 3 R> 
and feeding back ([37]) into ([T7]) we get 

EL - -o— 

whose integration with respect to r yields 



U = 



R 2 • 



where £ is an integration function of t, implying 
_ 3C(t) 



R 3 



(37) 



(38) 



(39) 



(40) 



In the case when the fluid fills all the sphere, including 
the centre (R(t, 0) = 0), we have to impose the regularity 
condition £ = 0, implying U = 0. Therefore if we want 
the expansion-free condition applied to all fluid elements 
to be compatible with a time dependent situation (U ^ 
0) we must assume that either 

• the fluid has no symmetry centre, or 



• the centre is surrounded by 
section of another spacetime, 
the rest of the fluid. 



a compact spherical 
suitably matched to 



Here we discard the first possibility since we are partic- 
ularly interested in describing localized objects without 
the unusual topology of a spherical fluid without a centre. 
Also, within the context of the second alternative we have 
chosen an inner vacuum Minkowski spherical vacuole. 

Thus the kinematical condition = is consistent 
with an evolving spherically symmetric fluid if there is a 
vacuum cavity surrounding its centre. 

We can also arrive at this latter conclusion by the fol- 
lowing qualitative argument. Since the expansion scalar 
describes the rate of change of small volumes of the fluid, 
it is intuitively clear that in the case of an overall expan- 
sion (contraction) , the increase (decrease) in volume due 
to the increasing (decreasing) area of the external bound- 
ary surface must be compensated with the increase (de- 
crease) of the area of the internal boundary surface (de- 
limiting the cavity) in order to keep vanishing. 
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IV. THE PERTURB ATIVE SCHEME 

We shall now describe in some detail the perturbative 
scheme which will provide us with the main equation 
required for our (in)stability analysis. 

We assume that initially the fluid is in static equilib- 
rium, which means that the fluid is described by quanti- 
ties only with radial coordinate dependence. These quan- 
tities are denoted by a subscript zero. We further sup- 
pose as usual, that the metric functions A(t, r), B(t,r) 
and R{t, r) have the same time dependence in their per- 
turbations. Therefore we consider the metric functions 
and material functions given by 



A(t,r) = 


Mr) + 


eT(t)a(r), 


(41) 


B{t,r) = 


-- B (r) + 


eT(t)b(r), 


(42) 


R(t,r) = 


-- Ro(r) + 


eT(t)c(r), 


(43) 


H(t,r 


= M r ) 


+ £}i{t,r), 


(44) 


P r (t,r) = 


-- Pro(r) - 


VtP r (t,r), 


(45) 


P±(t,r) = 


P±o{r) 4 


eP x (t,r), 


(46) 


m(t, r) 


= mo(r) 


+- em(t, r), 


(47) 




0(t,r) 


= eQ(t,r), 


(48) 




a(t,r) 


= £°(t,r), 


(49) 



where < e <C 1 and using the freedom allowed by the ra- 
dial coordinate we choose the Schwarzschild coordinates 
with Ro(r) — r. Considering (|41j|46p we have from (fT3l 
[TB)l for the static configuration 



1 



KflQ 



kP t o — 



n 



A" 

An 



(Bor) 2 
1 

A Q B 



2rf + BI-1 
A' 

a q 



A 



Bn 



(50) 
(51) 
(52) 



whereas from (|13M16[) we obtain for the perturbed quan- 
tities 



B 2 



\r) r\B ) \Bn r 



B' 



Bo 



Bo 
r 



b 

B~ 



ZkVoT—, 



(53) 



T 



A Bo 



(;)' 



b 



^0 



if 

TIT 



0, (54) 



T c 



T 



- 2 Af r +2 TB$ 



Ao 



b 

B~ 



2nP r nT 



Bo 



(55) 



T 



A 2 



b 

Bo 



kP±_ = To TT" + - + 



Ao 



T 



a 

~Aq 



b 

B~ 



-A 
Ao 



Bn 



2 4°-^ + - | ! — 



A 



B t 



(;)' 
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Ao 



2kP ±0 T 



Bo' 



(56) 



and for the expansion ([5]) and shear (fTTj) we have 



T 

Ao~ \Bo 



(58) 
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The Bianchi identities ([IS) and jUJ become with (|TO 
|46|) . for the static configuration 

P; o + (MO + P-o)^ + '(Pro - P±0) = 0, (59) 

A r 

and for the perturbed configuration, 



1 

A~o 



P + (mo + Pro)T— + 2(mo + P±o)T- 
B r 



0, (60) 



and 



1 

So" 



^ + ( M0 + Pr )T 



A 



(fi + Pr)f 



+2(P r0 -P i0 )T(^)' + 2(P r 



r 



0. (61) 



By substituting (j54"| into (|TO)) we can integrate it, and 
we find 



A* = 



6 c 
(Mo + fro)^- + 2(mo + Plo)- 
Po r 



P. 



(62) 



The total energy inside up to a radius r given by 
HI) with (|41M3j) and (gJJ becomes, 



T 



mo 
6 

Po" 



(64) 



Considering the second law of thermodynamics and the 
same arguments as given in , we can express a rela- 
tionship between P r and M given by 



P r = 



r0 



Mo + P-o 



M, 



(69) 



where Ti is the adiabatic index which measures the 
variation of pressure related to a given variation of 
density, thereby measuring the stiffness of the fluid. We 
consider it constant throughout the fluid distribution or, 
at least, throughout the region that we want to study. 
We recall that Ti coincides with the ratio of the specific 
heats for perfect Maxwell-Boltzman gas (30l - [32l |. 



V. NEWTONIAN AND POST NEWTONIAN 
TERMS 

Before dwelling to obtain the dynamical expansion- 
free equation, and since we are using relativistic units, 
the following comments should be helpful in order to 
identify the terms belonging to the Newtonian (N), post 
Newtonian (pN) and post post Newtonian (ppN) approx- 
imations. These terms are considered for the instability 
conditions stemming from the dynamical equation in the 
N and pN approximations. 

Thus, for N approximation we assume. 



From the matching condition (|30j) with (|45j) we have 



S (e) 



PrO = 0, P r = 



(65) 



For c ^ 0, which is the case that we want to study, with 
(|54|) . (|55|) and ([65]) we have 



T-aT = 0, 



(66) 



where 



Ao 
Po 



a 



4' 
An 



1 



b c 



(67) 



Solutions of (|66|) include functions which oscillate (cor- 
responding to stable systems) and those which do not 
(unstable ones). Since we are interested here in estab- 
lishing the range of instability, we confine our attention 
to the non oscillating ones, i. e. we assume that a(r), 
b(r) and c(r) are such that on r E ( e ), tt S ( C ) > 0. Then, 



T(t) = - exp (^(e) t) 



(68) 



representing a system that starts collapsing at t = — oo 
when T(— oo) = and the system is static, and goes 
collapsing, diminishing its areal radius, while t increases. 



Mo > -Pro, Mo > P±o- 



(70) 



For the metric coefficients, given in c.g.s. units, ex- 
panded up to pN approximation become 



An = 1- 



Gm 



Po = 1 



Gm 



(71) 



where Q is the gravitational constant and C is the speed 
of light. 

Next, from (|5ip with (|6"5|) we have (in relativistic units) 



A' KP r0 r 3 + 2m 
~An~ ~ 2r(r-2m ) ' 

and substituting into ([5§]) it yields 

KP, r or 3 + 2toq 



(72) 



P' 



(Mo 



p 



rO ) 



2r(r — 2mo) 
Writing ([73]) in c.g.s. units it becomes 



-(Pxo-Pro)- 
r 



(73) 



P' 



C-' 2 nPror a + 2m 



2r{r - 2C- 2 gm ) 



{H +C- 2 P r0 ) + -(P± -Pro). 

r 



Expanding (|74p up to terms of C 4 order we obtain 



(74) 
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Pro = -3 ^°™° + ^(P±o - Pro) - ^3 {2Q n ml + P r0 m r + ^fioP-ror 



4SVomo + 2gP r0 m 2 r + g^oPramar 4 + ^P r > 5 ) . (75) 



Hence from ([751) we have the following terms for the dif- 
ferent orders of approximations: 



N order: terms of order C°; 
pN order: terms of order C~ 2 ; 
ppN order: terms of order C~ 4 . 



(76) 
(77) 
(78) 



showing that the perturbed energy density of the system 
stems from the static background anisotropy. This fact 
is easily understood if we recall that under the expansion 
free condition, as it follows from (|20|) . changes in /i for any 
given fluid element, depend exclusively on the pressure 
anisotropy. 

On the other hand, with (l69l and (18211 we have 



VI. THE DYNAMICAL EXPANSION-FREE 
EQUATION 



P r = 2Ti P, '° (Pro - P±o)T-. (83) 
Mo + Pro r 



With the equations so far obtained we can build the dy- 
namical equation that we specialize to the expansion-free 
condition which is the aim of our study. The obtention 
of the dynamical equation is done via (|6ip . 

The expansion-free condition = implies from (|57p 



B r' 



(79) 



which together with the adiabatic condition (|54|) pro- 
duces, 



T ( r 2 c 



r 3 B \Ao 



0, 



(80) 



implying 



Ao 



• = fc 7T> ( 81 ) 



where k is a constant. With (l79l we have for J62 



Jx = 2(P r0 - P ±0 )T-, 
r 



Two further useful relations we can obtain. One from 
(ESI), 



An 



1 



A 



Ho + Pi 



,0 



PrO 



'-(PrO-P±o) 



and another from ([50|) and (|63|) 

" 3 — 2m 



B' 



B 2r(r ~ 2m ) 



(84) 



(85) 



We want to study the instability conditions for the 
expansion-free fluid up to the pN approximation, since in- 
termediate calculations are rather long we shall put them 
in an Appendix. 

Considering the N approximation of the dynamical 
equation (jM|) by using (llOOp and (I104j|106p and that 



terms with P r o / Ho being of ppN order it reduces to 

3^0 + n\P r0 \r + 2 - ^ ^ = 4 K (5P r0 - 2P ±0 ), 

(86) 

(82) or, using(|24| and rearranging terms 



18' 



IP' 



rO 



4 , a £<=> 2 

r H uiov 



2 k 



(5P r0 - 2P ±0 )r 3 




r 2 dr 



fJ-or 



(87) 



where we assume P/ < 0- 



For the pN approximation we have 
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-6 a 



*£(<=) 



\ r / r 



Let us first consider the N approximation. The first ob- 
servation to be made is that in the absence of a single pa- 
rameter such as Ti the assessment of the instability range 
depends in a rather complicated way on different struc- 
tural properties of the fluid, such as pressure anisotropy 
and the radial profile of the energy density. 

Indeed, for the onset of instability we need (|87|) to be 
satisfied, and since the two terms at the left of (|87|) are 
positive, then instabilities may develop only if so is the 
combination of the two terms at the right of (|8"T)l . For 
that to happen it will be sufficient that P r o > (2/5)P±o 
and that the last term in ([87]) be positive. Let us explore 
this possibility in some detail. 

Thus, let us consider an energy density profile of the 
form /io = /3r™ where (3 is a positive constant, and n is 
also a constant whose value ranges in the interval — oo < 
n < oo. In this case the last term in ([57)1 becomes (for 



K0 



„"+3 



6(3 + n) 
whereas for the n = 

IT 



4-n-7 



3+n 



-3 case we obtain 



7 log 



?"£(>) 



1 



(89) 



(90) 



Then the following possibilities arise: 
1. n < 0, n ^ -3. 

In this case (|8"9"|) will be positive if 

l/(n+3) 



r > TV 



(91) 



Thus the maximal range of instability decreases 
from its value for |n| close to 3, for which we have 
(see case 3 below) 



to 



r > r E(i) e 



r > r S (i) 



1/7 



1/3 



r E( i)1.15 



r sW 1.20 



(92) 



(93) 



corresponding to n = (the incompressible fluid). 
On the other hand as \n\ increases from |n| =3 
the maximal region of instability will embrace the 
whole fluid, since 



lim 

n— > — oo 



7 



4 - n 



l/(3+7l) 



= 1 



(94) 



In other words, for n < 0, the range of instability 
increases as the absolute value of n increases. 

2. n > 0. 

In this case the positivity of the last term in (|87|) 
is, again, assured in the region 

i/(n+3) 



r > r EW 



7 



4 — n 



(95) 



The range of instability decreases with n, vanishing 
for n > 4. 

3. n = —3. 

In this particular case, the positivity of the the last 
term in (|87[) is assured in the region 



r > r E(<) e 



1/7 



r s(l ,1.15 



(96) 



Of course, if P r o < (2/5)Pj_o (for the same energy den- 
sity profile), the maximal region of instability diminishes. 

At pN order the situation is essentially the same, with 
the relativistic effects at first order taken into account. 

Observe that Ti does not play any role in the N and pN 
orders. In other words at N order, no matter how stiff 
is the material, the system will be unstable as long as 
is satisfied. This is at variance with the result in the 
non- vanishing expansion case when it appears that the 
range of instability is defined by Ti < (4/3)+ anisotropic 
term (see Q for details). A similar remark applies for 
the pN order. 

As mentioned before the fact that Ti does not enter 
into and (|8"8"|) becomes intelligible when we recall 
that the expansion-free collapse (close to the Newtonian 
regime) proceeds without "compression" of the fluid. Ac- 
cordingly the stiffness of matter is irrelevant for the on- 
set of instabilities, the later being dependent on the local 
anisotropy of pressure and energy density inhomogeneity. 

It should be observed that in the process of collapse 
(under the expansion-free condition) the energy density 
may increase or decrease depending on the difference be- 
tween the velocity of the inner and the outer boundary 
surface. This diference in turn is determined by the local 
anisotropy of pressure and the radial profile of the energy 
density. 



VII. CONCLUSIONS 

We have seen so far that in the study of dynamical 
instability, under the expansion-free condition, at, both, 
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the Newtonian and post Newtonian regimes, the range of 
instability is defined by the local anisotropy of pressure 
and the energy density radial profile, but not by the adia- 
batic index I\. This implies that in the above mentioned 
regimes the stiffness of the fluid, measured by Ti, which 
generally plays the central role in the definition of the 
instability range, is irrelevant here. This fact strength- 
ens further the relevance of local anisotropy of pressure 
and energy density inhomogeneity in the structure and 
evolution of self-gravitating objects. 

It should be stressed that any possible model is further 
constrained by physical requirements such as positivity 
of energy density, sound speed less than light (energy 
density greater than pressure) , and stability of local os- 
cillations modes. 

The association of a cavity with the expansion-free 
fluid distribution implies that the above presented study 



describes the instability range of the cavity (keeping the 
expansion- free condition). We observe that the pertur- 
bation of the cavity itself produces a similar equation to 
T(t) as in (1551) and it imposes a further constraint to the 
system a s(e) = a E <« . 

It is worth mentioning that the role played by the 
anisotropy and the energy density inhomogeneity, as it 
follows from the previous section, is fully consistent with 
the results obtained from the study on the influence of 
those two factors on the active gravitational (Tolman) 
mass [II], presented in 34] -[36]. 



Indeed, for the Tolman mass mj- [33| interior to a 
sphere of radius r, up to pN approximation, the following 
expression may be obtained in the static or quasi-static 
case (see equation (47) in 34], or equation (58) in [35[ or 
equation (50) in [36| and the discussion therein) 



J 



m T = (m T ) s(c 



/ (PrO-Pxo)--^ / ^df 



dr. 



(97) 



r 



where (fTTj) has been used. 

From the above it is evident that fi' < and P r o > 
P±o increase the Tolman mass. That result, together 
with the fact that the Tolman mass may be interpreted 
as the active gravitational mass, provide full support to 
the conclusions on the instability range obtained in the 
previous section. 

It is worth noticing that if we keep the expansion-free 
condition all along the collapse, then as soon as the fluid 
approaches the centre there will be a blowup of the shear 
scalar, as implied by (|40l) . But it is shown in [l5j that 
sufficiently strong shearing effects near the singularity 
delay the formation of the apparent horizon, implying 
the appearance of a naked singularity. In other words 
the expansion-free condition provides a simple scenario 



for naked singularity formation. The obvious relevance 
of such an effect strengthens further the interest of the 
problem discussed here. 

Finally, let us mention that an extensions of these re- 
sults for f(r) gravity theory, have been recently presented 

IB, M. 



Appendix 

Inspecting (|53")) we see that P r is of ppN order, as well 
as p,A' /A by considering ([8"2]l and (|53")l . hence (1BT1) to- 
gether with and (|56"j) reduces to 



J 



K(fi + P r0 )r [ —J + 2n(P r0 - P±o)r ^ 



8kP ±0 - 



a 

a7> 



(;)' 



A'o B' , 1 



A B r 



(;)' 



a 



A B r 

_ p Q^) c 
Al r 



0. (98) 



r 



From flSSJl with ([6"6"j). ([79"| and ((5TJ we have 



A' 



2 A' 



lr R 2_^ fB°\ 2 ] 



(QQ\ 



The first three terms of with flU) and up to pN 
order become 
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k(^o + P r o)r 
+ (^o + P r o)r 



a 

Ao~ 



, / c \ ' c k/c An 

2«(P r0 - Pxo)r - - 8 K P ± o- = 

\r/ r r 



Un^ProrBl + 2P.;, Q + ^(5P r0 - P±) 



+ (^o + P r o)r 



nk 

= J { 2K^ P r or + [ 1 



,m 



a*., (1 + 3^) 



_ m ° '\ 
r ) 



2P r ' +-(5P r0 -P±o) 



6 15 
r"2 ~ Y 



(100) 



Now we calculate the following terms, 



a 



A' 



B' 



A B 



a 



A 



1_A 

r 2 An 



£ (5-9B 2 ) 



Bo 



(*o 2 + l) + >o 2 -l) 



-2kP t0 Bq 



Bo 
Ao 



4' 
A 

An 



_ 2Aj,Pj, 
r Ao B 
B^ l 
Bo r 



(101) 



and 



(;)' 



A'o B' 4 
A B r 



— k ^ 



4. 

A 



Ao^o 



11 A' 



3 B 



A Bo r A r B 



(102) 



where (JHTJ and (|99|) have been used. With (fTUTj) and 
(fTU2]) we can build 



R2 

-°o 



a 

Ao" 



2 A P + rJ UoJ V Ao B + r) \r) 



= Ank^- 



PrO 



- 6k 



Ao 



1 A 



( 3^ + 2)-^--B B + -(B 2 -l) 



B5V 3 



Aj, 
A 



A^Po 
A B 



0%^ 


fA 




A r 2 


Uo 





(103) 



By using (fTTj) , (|84|) and ff85|) we obtain the following ex- 
pressions up to pN order, 



Ank^ 



PrO 



Bo 



Auk 



Pn 



= ( p; - ^ + >/, Po - ^p;, 



(104) 
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-6k 



B V 



A 



+ 5 



4' R' 
A Bo 

r / no 

9^1 — 

r J r^o 



(SBt + ^-lBoB'+^Bl-l) 



Kr 3 /j, 



Mo 
Mo 



2mo 



-(PrO-Pxo) 



2 



2r 2 M o 

2' 



PrO + _ (-PrO - P±o) 

r 



1 m 3 /m \ 
r 2 V r / 



+ 



5 (v)' 



,(105) 



-2 



^or 2 Uo 



5 \ 
p ; 




{- 


K 

+ 2 


1 + 3^ + 
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I 



mo\ J_ 
r ' Mo L 

2" 



Pj-O 



-(PrO-Pxo) 

r 



v r y 



1 

»"Mo 

r 



mo +3 /moV 
r \ r J 



(106) 



some terms of ppN order appearing in some of the equa- VI. 
tions above will be excluded in the analysis of section 
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